RD-A146  686  THE  STflBILITV  OF  AN  INFINITE  ELASTIC  PIPE  WITH  UNIFORM  i/1 
MEAN  FLOW (U)  ADMIRALTV  MARINE  TECHNOLOGV  ESTABLISHMENT 
TEDDINGTON  (ENGLAND)  E  A  SKELTON  FEB  84 
UNCLASSIFIED  AMTE(N)/TM8405?  DRIC-BR-92675  F/G  20/4  NL 


*->t  v  y  l.  u  r  J 


UNLIMITED 


AKFE(  H  )TNB40S7 


THE  STABILITY  OP  Ml  INFINITE  ELASTIC 
PIPE  WITH  UNIFORM  MEAN  PLOW 


BY 


E  A  SKELTON 


The  pipe  conveying  Incompressible  fluid  Is  node lied  using  linear  bean 
thfcory.  The  stability  of  the  systee  Is  analysed  using  a  causal  approach  which 
is  able  to  distinguish  between  absolute  and  convective  Instabilities.  Closed- 
fora  expressions  for  the  stability  boundaries  of  the  bean  theory  aodel  are 
derived  analytically.  The  dissipation-free  pipe  Is  absolutely  unstable  if  the 
ratio  of  the  aass  per  unit  length  of  the  pipe  to  that  of  the  contained  fluid 
is  greater  than  one-eighth,  otherwise  it  is  convectlvely  unstable  at  low 
excitation  frequencies  and  stable  at  high  frequencies.  The  pipe  is  found  to 
be  absolutely  unstable  when  doping  is  included.  Some  numerical  results  are 
presented  for  a  steel  pipe  and  a  rubber  pipe,  both  of  Which  contain  water 


ANTE  (Teddington) 
Queen's  Road 

TEDOIMGROM  Kiddx  TW11  OIN 


42  pages  February  1984 

12  figures 


© 

Copyright 

Controller  HMSO  London 


UNLIMITED 


COHTEHTS 

1.  Introduction 

2.  Beam  Theory 

3.  Mathematical  Analysis  of  Beam  Stability 

(a)  General 

(b)  Coalescing  poles 

(c)  Response  to  Time-harmonic  Switch-on 

(d)  Effect  of  Damping 

4.  Numerical  Results 

(a)  General 

(b)  Beam  Theory  Results 

5.  Concluding  Remarks 

Figures  1-12 


Appendix  Causal  Approach  to  the  Stability  Analysis  of  Linear  Systei 


LIST  OP  SYMBOLS 


Wr(S.t) 

Pr(z.t) 

W(*,t) 

P(*.t) 


W(  z,u) 


W(  «, « ) 


D<a,u) 


®»h,ps 


Instantaneous  value  of  beam  displacement 


Instantaneous  value  of  bean  excitation 


displacement  amplitude 
excitation  amplitude 
Tapi ace  transform  of  W(z,t) 
Fourier  transform  of  W( z,  w) 


dispersion  relation 


axial  and  circumferential  wavenumbers 


frequency  in  Hz 

radian  frequency  (-2tt£) 

denotes  the  line  of  Laplace  integration 

density,  sound  velocity  and  axial  flow  velocity  of 
interior  fluid 

Mach  number  of  fluid,  Mj-Oj/Cj 
radius,  thickness  and  density  of  shell 
Young's  modulus  of  the  shell 


Poisson's  ratio  of  the  shell 


hysteretic  loss  factor  of  the  Shell 
mass  per  unit  length  of  fluid,  M-ira2pj 
mass  per  unit  length  of  beam,  m»2trahps 
area  moment  of  inertia  of  beam,  I«wha3 


(as  a  superscript)  denotes  complex  conjugate 


mraoDocTZON 


Much  theoretical  work  [1-9]  has  been  dona  on  the  vibration  and  sound 
radiation  of  fluid- loaded  structures ,  such  as  plates  and  shells,  excited  by 
tine— harmonic  — chanlcal  forces  or  acoustic  point  sources.  These,  and  other, 
authors  deal  with  wave  propagation  In,  and  vibration  response  of,  systems  In 
tdiiCh  the  fluid  Is  neeimmd  to  be  stationary.  Practical  pipework  systems, 
however,  usually  convey  fluid  whose  action  aay  significantly  Change  the 
dynamic  characteristics  of  the  systea.  In  particular,  the  interchange  of 
energy  between  the  moving  fluid  and  the  pipe  aay  result  In  an  Instability 
which  causes  excessive  vibration  and  sound  radiation. 

The  stability  analysis  of  systems  conveying  fluid  usually  proceeds  via 
the  system  dispersion  relation.  Complex  values  of  frequency  corresponding 
to  real  values  of  wavenumber  are  obtained  and  the  sign  of  the  imaginary  part 
of  these  frequencies  indicates  whether  or  not  the  system  Is  unstable. 

Dowell  [10]  has  used  this  approach  to  investigate  the  flutter  behaviour  of 
Infinite  shells.  Be  concludes  that  the  shall  theory  used  was  inadequate  for 
the  n-1  harmonic,  but  for  each  n*l  circumferential  harmonic  a  flutter 
instability  occurs  above  a  critical  flow  speed,  the  lowest  such  speed  being 
the  one  associated  with  the  n*2  harmonic.  The  effect  of  damping  In  the  shell 
was  not  considered. 

Mslcher  [11]  describes  an  alternative  approach  to  stability  analyses 
which  was  used  by  Briggs  [12]  In  connection  with  beam-plasma  interactions. 

This  approach  demands  that  the  system  Is  causal,  that  is,  it  cannot  respond 
until  an  excitation  Is  applied.  At  a  sufficiently  long  time  after  the 
excitation  has  been  applied  three  types  of  behaviour  nay  be  identified. 
Firstly,  there  Is  the  stable  configuration  in  which  the  response  is  finite 
everywhere  In  space-time.  Secondly,  there  is  the  correctively  unstable 
configuration  in  which  the  response  increases  exponentially  with  distance  from 
the  excitation.  Thirdly  there  is  the  absolutely  unstable  configuration  In 
which  the  response  at  all  points  In  the  system  Increases  exponentially  with 
time.  Brazier-Smith  a  scott  [13]  have  applied  this  causal  approach  to  the 
problem  of  uniform  Incompressible  flow  over  an  Infinite  undamped  elastic 
plate.  The  excitation  was  assumed  to  be  a  time  harmonic  force  switched  on  at 
t-o.  They  found  that  at  low  excitation  frequencies  there  are  convective 
Instabilities,  and  that  above  a  critical  flow  velocity  there  axe  also  absolute 
instabilities.  Atkins  [14]  has  confined  the  main  results  of  this  paper  and 
has  also  demonstrated  that  plate  damping  may  destabilise  some  free  waves  at 
all  flew  speeds  and  does  not  remove  the  absolute  instability  which  occurs 
above  the  critical  velocity. 

In  the  present  work  the  causal  approach  to  stability  determination  is 
applied  to  an  infinite  pipe  with  uniform  axial  fluid  flew.  The  pipe  la 
modelled  as  a  beam  conveying  an  Incompressible  fluid.  The  effects  of  pipe 
damping  are  Included  in  the  model.  The  stability  Characteristics  of  the  ben 
theory  model  are  deduced  analytically,  whereas  a  numerical  approach  would  be 
necessary  to  determine  the  stability  Characteristics  of  shell  theory  models. 


2.  BEAM  THEORY 


TIM  transverse  displacement  wr(  z,t )  of  a  pipe  conveying  an  inviscid 
incompressible  fluid  is  often  assumed  to  satisfy  the  differential  equation 
of  beam  theory  [15],  viz., 

zia\ ^/dx4  +  Mo**2* ^az2  +  2wla2wr/azat  +  (M*m)a2wr/at2  -  Pr(z,t)  (2.1) 

Where  M  is  the  mass  per  unit  length  of  the  fluid,  m  is  the  mass  per  unit 
length  of  the  pipe  and  I  is  the  area  moment  of  inertia  of  the  pipe. 

The  causal  solution  for  the  displacement  may  be  obtained  by  expressing  it 
as  the  Laplace-Pourier  transform 

otfier  oo 

W(z,t)  -  (V4ff2)  |  j  | «*P( iaz-iut )dodu  (2.2) 

-ootiff  -« 

Where  the  dispersion  relation,  D(a,u»),  is  obtained  from  equation  (2.1)  as 

D(a,u)  -  Ela4  -  IBJ2o2  +  2*0^  -  (Mtm)u2  (2.3) 


3.  MATHEMATICAL  ANALYSIS  OP  BEAM  STABILITY 


(a)  General 


The  symmetrical  nature  of  the  dispersion  relation,  viz., 

0(a,u)  «  D*(-a*,-«  )  (3.1) 

means  that  it  is  necessary  to  consider  only  those  values  of  u  such  that 
Re[u]»0,  the  stability  analysis  with  Re[u]<0  being  almost  identicals  the 
root  loci  in  the  complex  a-plane  obtained  for  those  u  with  Re[u>]<0 
are  reflections  in  the  imaginary  a-axia  of  those  obtained  for  Rs[<i>]>0. 

The  analysis  of  beam  stability  is  further  simplified  by  the  absence  of 
branch  cuts  in  the  complex  a-plane. 

The  first  step  in  the  stability  analysis,  as  described  in  the  Appendix, 
is  to  search  for  coalescing  poles  of  the  dispersion  relation  which  originate 
in  opposite  half-planes.  If  such  poles  exist,  they  indicate  the  presence  of 
absolute  instabilities  for  both  Impulse  and  time  harmonic  switch-on 
disturbances,  secondly,  in  the  absence  of  coalescing  poles  the  form  of  the 
response  to  the  time-harmonic  disturbance  is  determined,  the  response  to  the 
impulse  being  bounded  in  both  space  and  tlma. 


(b)  Oo*l— dug  Pol— 


Coal— icing  polos  aro  px—ont  at  tho—  values  of  Up  and  otp  which 
satisfy  the  simultaneous  equations 

D(a,u>-0  (3.2) 

9D(a,u)/dcrt)  (3.3) 

It  is  straight forward  to  solvs  tbs—  aquations  fox  tbs  beam  to  giva  tba  two 
solutions 

«p-oJ(  2M/EI)1/2(  l*v(  1— BSl/tl)  )  1/2(  3*V(  l-8m/M)~3/2  (3.4) 

o-O  (H/2EI)1/2(  1*V(  1-es/M) )1/2( 3*v( l-8a/M)_1/2  (3.5) 

p  1 

whara  fee  tba  first  solution  tba  positive  signs  are  taken  throughout,  and 
for  tba  second  solution  tbs  negative  signs  axe  taken  throughout,  lbs 
r—nlnlng  two  solutions  of  equations  (3.2)  and  (3.3)  axe  of  opposite  sign 
to  those  given  above,  and  in  particular  have  Re[<dp]<0.  For  causes  where 
8n/N4l  both  of  the  solutions  for  Up  given  by  equation  (3.4)  are  purely 
real,  and  therefore  represent  transient  solutions,  even  If  the  root  loci  of 
the  coalescing  poles  originated  In  opposite  half-planes. 

If,  however,  8m/M>l,  then  the  — oond  solution  for  (Up,ap)  is  such 
that  up  has  a  positive  imaginary  part,  indicating  that  an  absolute 
Instability  is  present  if  the  root  loci  of  the  coalescing  poles  originated 
in  opposite  half-planes .  me  corresponding  value  of  <Xp  has  a  positive 
real  part  and  a  negative  imaginary  part.  She  presence  of  an  abeolute 
Instability  is  determined  —  fellows. 

(I)  At  any  real  frequency  Uj,  the  four  root  loci  for  u-o^+io,  where  a 
iM  large  and  positive,  are  found  fxosi  the  various  caahinations  of 

a  »  [o2(IH«)/4EI)1/4(*l±i)  (3.6) 

Hence,  each  of  the  root  loci  originates  in  a  separate  quadrant  of  the 
coupler  e-plane. 

(II)  In  order  to  determine  whether  ox  not  the  coalescing  pol—  originated 
in  the  same  half-plane,  it  is  convenient  to  solve  the  dispersion  relation, 
equation  (3.2),  for  u,  viz., 

«-a[HDl*(EI(bf«)a2-— «J2) 1/2  ]/(»*♦«)  (3.7) 

and  to  co— ider  where  the  root  loci  nay  cross  the  real  and  imaginary  a-aa— . 

(ill)  If  a  root  locus  crosses  the  imaginary  a-axis,  *:  a-ir  say,  then 
from  aquation  (3.7) 


u-lr[IK;i*i(Bl(  N4«  )r24«K;2 } 1/2  ]/(  MHa ) 


(3.8) 


and  because,  by  definition,  on  &  root  locus  lm[w]  vast  bo  positive, 
equation  (3.9)  shows  that  r  Must  also  be  positive.  Hence,  no  root  locus 
secy  cross  the  negative  imaginary  a-axis. 

(iv)  If  a  root  locus  crosses  the  real  a-axis,  at  «  say,  then  by 
considering  the  real  and  iaaginary  parts  of  equation  (3.7) 

(Mfa)EXs2~fl*cr^<0 

and 

»r-sw /(*♦«) 

Since  the  analysis  has  been  restricted  to  only  those  values  of  u  with 
positive  real  part,  equation  (3.10)  shows  that  s  oust  also  be  positive. 
Hence,  no  root  locus  may  cross  the  negative  real  a-axis.  nils  result, 
together  with  the  result  of  paragraph  ( iii )  shows  that  the  locus  which 
originates  in  the  lower  half-plane,  with  Re[a]<0,  oust  resain  in  that 
quadrant  of  the  coaplex  a-plane,  and  in  particular  it  say  not  coalesce 
with  any  other  loci,  thus,  if  the  coalescing  poles  originate  in  the  sane 
half-plane,  it  suet  be  the  upper  half-plane. 

(v)  Suppose  that  the  loci  of  the  coalescing  poles  both  originate  in  the  upper 
half-plane,  then,  because  the  iaaginary  part  of  ctp  is  negative,  they  oust 
both  cross  the  (positive)  real  cx-axis.  Equations  (3.9)  and  (3.10)  show 
that  for  a  given  real  frequency  c^,  loci  nay  cross  the  real  ct-axis  at 

(at  sost)  one  point.  By  substituting  this  (real)  value  of  a  into 
equation  (3.7)  it  is  clear  that  only  one  (positive)  value  of  a  is  possible 
for  both  loci  at  this  point.  Hence,  coalescence  occurs  on  the  real  a-axis 
and  a  contradiction  to  the  results  of  equation  (3.5)  has  been  obtained. 

Thus  the  coalescing  poles  nay  not  both  originate  in  the  upper  half-plane. 

Together,  the  results  of  paragraphs  (iv)  and  (v)  shew  that  the  loci 
of  the  coalescing  poles  originate  in  opposite  half-planes,  and  therefore 
indicate  that  the  pipe  is  absolutely  unstable  if  8m/M>l,  and,  as  described 
in  the  Appendix,  the  pipe  response  is  dominated  by  the  frequency  Up. 


(c)  Response  to  Time-harmonic  Switch-on 


in  the  absence  of  an  absolute  instability,  the  fore  of  the  pipe  response 
to  a  forcing  excitation  at  a  particular  frequency  uq  soy  be  determined.  Fo: 
8*yHCl  there  is  no  absolute  instability  present,  and  the  behaviour  of  the 
pipe  depends  on  the  forcing  frequency.  Zt  follows  from  equations  (3.9)  and 
(3.10)  that  a  locus  for  u-ttQ+io'  crosses  the  real  a-axis  only  if 

«0<H02(»lfmrl[iWEI(  *♦«)] 1/2  (3.11 

this  is  a  necessary,  but  not  sufficient,  condition  for  a  convective 
instability  since  the  locus  say  subsequently  return  to  the  real  a-axis, 
in  trtiicft  case  the  pole  would  represent  a  propagating  wave. 

By  considering  the  dispersion  relation,  in  the  form  given  by  equation 
(3.7)  for  real  values  of  a  and  u  it  is  obvious  that  the  branch 


(3.9) 

(3.10) 


(3.12) 


w(a)"a[Mtlj+(EI(IHm)a  -sm^}  '  J/(H4«) 

is  a  wonotonic  incrsasing  function  of  a.  in  ths  regions  of  the  real  a-axis 
where  it  exists  as  a  zsal  function.  Clearly,  however,  the  branch  given  by 

«K  a  V-aCW^tEKUHi^-sHtJ*} 1/2  ]/(*+«)  (3.13) 

is  not,  and  the  equation 


u0«u(a)  (3.14) 

■ay  have  2,  3,  or  4  real  roots  a,  depending  on  the  value  of  uq.  The 
boundaries  between  these  intervals  on  the  w-axia  are  determined  by  solving 
the  equation 


du/0o"O  (3.15) 

Where  u(a)  is  given  by  equation  (3.13).  The  pipe  behaviour  may  then  be 
descibed  in  each  of  three  frequency  regions,  using  arguments  similar  to 
those  used  in  section  3(b).  The  behaviour  is  as  follows! 


Por 


«0<O*(  2M/E1 ) l/2( 1— V( 1-8S/M)  ) l/2( 3— V( l-9m/H)~3/2  (3.16) 


the  pipe  is  convectively  unstable  in  the  region  s>0.  in  this  region  two 
waves  Whose  phase  propagates  downstream  from  the  excitation  exist,  one  of 
which  decays  exponentially  with  distance,  whilst  the  other  grows  exponentially 
with  distance.  In  the  upstream  region,  z<0,  there  are  two  free  waves  which 
propagate  without  decay,  one  whose  phase  propagates  downstream  and  the  other 
upstream.  However,  the  group  velocity,  9u/0a,  of  both  of  these  waves  is 
negative,  indicating  that  in  both  cases  the  energy  flows  upstream  from  the 
source  of  excitation. 


Por  “o>0i(  2>VEI)1/2(  W*H) )  l/2(  3— v(  l-8m/*)-3/2 


(3.17) 


and 


«X(  ZM/E1 )  l/2(  l+v(  1-8W/M )  ) l/2(  3-M  l-9m/H)~3/Z 


the  pipe  behaviour  is  stable  and  consists  of  four  free-waves.  In  the 
upstream  region,  s<0,  the  phase  of  one  wave  propagates  upstream  from  the  point 
of  excitation,  whilst  ths  phase  of  the  other  propagates  downstream  from 
infinity.  Again,  in  both  cases,  the  group  velocity  is  negative  and  the 
energy  flow  is  upstream  from  the  point  of  excitation,  in  the  region,  s>0,  the 
phase  of  both  free  waves  propagates  downstream  from  the  source.  The  group 
velocity  of  both  of  these  waves  is  positive,  and  the  energy  flow  of  each 
wave  is  downstream  from  the  source  of  excitation  unless 

wft<  MtJ2(i»fm)"l(iea/si(ilH»>]1/2 


when  the  group  velocity  of  on*  of  the  «ma  is  negative,  and  the  energy  flow 
associated  with  that  wave  Is  upstream  from  infinity. 

For 

»0>0*i2M/El)l/2(l+V(l-am/H))i/\3+V(l-9m/H)-3/2  (3.18) 


tha  plpa  behaviour  is  stable  and  consists  only  of  propagating  and  decaying 
waves.  In  the  downstream  region,  s>0,  two  waves  travel  downstream  from  the 
source,  one  of  Which  is  a  free  wave  which  propagates  without  decaying  and 
the  other  is  evanescent  and  decays  exponentially  with  distance.  In  the 
upstream  region,  z<0,  one  free  wave  propagates  upstream  without  decaying. 

The  amplitude  of  the  second  wave  in  the  z<0  region  decays  exponentially  with 
distance  from  the  excitation  and  is  therefore  stable. 


(d)  Effect  of  Damping 


A  positive  hysteretic  loss  factor,  q,  may  be  introduced  Into  the 
mathematical  model  by  making  E  complex,  vis.. 


E-E<l-in). 

E-*(H-iT)), 


for  w>0 
for  u<0 


(3.19) 


where  i)  is  small  compared  with  unity.  The  symmetry  of  the  dispersion 
relation,  equation  (3.1),  is  preserved  which  means  that  it  is  again  sufficient 
to  consider  only  those  values  of  u  with  Ra[u]>0. 

Equations  (3.4)  and  (3.5)  again  give  the  position  of  coalescing  poles 
except  that  now,  for  small  n. 


«  mu  (1+1V2) 
P  P 

a  «  (  I+Itt/2  ) 
P  P 


(3.20) 


For  the  case  aa/M>i,  the  addition  of  damping  has  altered  the  location 
of  the  coalescing  poles  slightly,  but  not  enough  to  stabilise  tha  system 
from  its  previously  absolutely  unstable  state. 


For  the  case  S^IKl,  where  there  was  no  absolute  Instability,  the 
addition  of  damping  alters  tha  position  of  the  root  loci  slightly,  and  may 
result  In  coalescing  poles.  In  the  absence  of  damping,  both  solutions  for 
the  frequency,  given  by  equation  (3.4),  at  which  poles  coalesce  were  real  and 
thus  represented  transient  solutions.  Equation  (3.20)  Shows  that  with  damping 
both  of  these  solutions  have  a  positive  Imaginary  pert,  and  may  therefore 
represent  absolute  instabilities  if  the  coalescing  poles  originated  In 
opposite  half-spaces.  The  behaviour  of  the  root  loci  in  the  undamped  case 
•how  that  at  the  lower  frequency  the  coalescence  is  between  loci  which  both 
originate  In  the  upper  half  of  the  complex  a-plane,  but  at  the  higher 


Up-ty  2M/Z1 )  '  *(  l+v(  1-Bm/M)  )  '  (  3+V(  l-8m/H)  ) 


0-21) 


the  coaleeoence  is  between  one  pole  whose  locus  originates  in  the  lower  half¬ 
space  and  one  whose  locus  originates  in  the  upper  half-space.  There  is  thus 
an  absolute  instability  at  the  frequency  Wp  given  by  equation  (3.21). 
However,  because  the  exponential  growth  rate  is  proportional  to  n»  which  is 
small,  this  contribution  to  the  displacement  may  not  dominate  the  solution 
until  t  becomes  very  large.  For  small  or  moderate  values  of  the  tine  t, 
the  response  will  appear  to  be  approximately  as  was  determined  for  the 
undamped  case. 


4.  NUMERICAL  RESULTS 


(a)  General 


Fortran  programs  have  been  written  for  calculating  the  complex  roots, 
a,  of  the  dispersion  relation  given  by  equation  (2.3)  for  complex  values  of 
the  frequency  parameter  teei^+ic.  The  programs  are  written  in  double 
precision  complex  arithmetic,  which  is  simulated  [16]  by  the  use  of  double 
precision  real  arrays  of  leading  dimension  two.  This  was  necessary  because 
the  programs  were  run  on  a  PDP-11  computer  whose  single  precision  arithmetic 
word  length  of  32  bits  is  inadequate  for  this  type  of  problem. 

The  material  and  geometric  constants,  in  SI  units,  which  were  used  in 
the  calculations  for  Figures  6-11  are  as  follows i 

Steel  Pipe!  E- 19. 5x10 10  v-0.29  p-7700.0  a-l.OxlO-2  h-O . 05x10“ 2 


Rubber  Plpet  E- 3.23x10 10  v-0.4  p-1100.0  a-l.OxlO-2  h-O . 05x10" 2 


Water:  p- 1000.0  c- 1500.0  0-10.0 

Damping  in  the  pipe  wall  was  included  in  the  calculations  for  Figures  10  and 
11  by  setting  E,  the  Young's  modulus,  to  the  complex  value  E(l-in),  where 
the  numerical  value  of  n,  the  hysteretic  loss  factor,  was  chosen  as  0.2. 

This  value.  Which  is  an  order  of  magnitude  higher  than  practical  levels,  helps 
to  clarify  the  effect  of  damping  without  significant  Change  to  the  physical 
interpretation . 


(b)  Beam  Theory  Results 


Figures  6a  and  6b  Show  the  real  and  Imaginary  parts,  respectively, 
of  the  wavenumber-frequency  plot  for  the  undamped  steel  pipe,  obtained  from 
the  dispersion  relation  for  a  beam,  equation  (2.3).  The  branches  labelled 
1  and  2  are  purely  real,  and  branches  3  and  4  are  complex  conjugate  roots. 


The  positive  real  branch,  labelled  1,  has  a  mini—  vaJ.ua  at  approximately 
0.198s  Which  indicates  an  infinite  group  velocity  and  also  a  possible 
instability  at  this  frequency.  Clearly,  because  these  brandies  of  the 
dispersion  relation  do  not  intersect,  no  poles  may  coalesce  at  real  values 
of  the  frequency. 

Figure  7  Shows  a  root  locus  stability  plot  for  the  four  roots  (labelled 
a,  b,  c  and  d)  of  the  undamped  steel  pipe,  all  the  loci  labelled  1  correspond 
to  the  frequency  WJ./2TT-0.04H*.  The  loci  labelled  la  and  lb  terminate  mi  the 
real  cr-axis  and  for  a  time-harmonic  switch-on  excitation  at  this  frequency 
therefore  represent  free  waves.  Whose  phase  propagates  upstream  from  the 
excitation  and  downstream  from  infinity,  respectively.  The  energy  of  both  of 
these  waves  propagates  upstream  from  the  excitation.  The  locus  labelled  lc 
crosses  the  real  a-aads  and  terminates  in  the  lower  half-plane.  This 
represents  a  convective ly  unstable  wave  propagating  dcwnatream  from  the 
excitation.  The  locus  labelled  id  terminates  in  the  upper  half-plane,  and, 
because  there  are  no  branch  cuts  in  the  complex  a-pl&ne,  may  be  identified 
with  an  evanescent  wave  propagating  downstream  from  the  excitation. 

♦ 

The  loci  labelled  2  and  3  correspond  to  frequencies  Up^tr-o. 08H*  and 
<ar/2ir-o.l2Hs,  respectively,  and  represent  responses  similar  to  those  at  0.04HX 
The  loci  labelled  4  and  5,  corresponding  to  frequencies  ur/2rr«0.16Hs  and 
ur/2fpK).208B,  however  show  very  different  behaviour.  The  locus  labelled  4b 
terminates  in  the  lower  half-plane,  and  may  be  identified  with  a  wave  whose 
phase  and  group  velocities  both  propagate  downstream  from  infinity.  The 
amplitude  of  this  wave  decays  exponentially  with  distance  from  the  excitation 
point,  s-O,  and,  according  to  NslCher's  classification  it  must  therefore  be 
described  as  stable,  in  fact  the  amplitude  of  this  wave  increases  in  the 
direction  of  propagation,  i.e.  it  appears  to  be  a  convectively  unstable  wave 
propagating  downstream  from  infinity.  The  locus  labelled  4c  terminates  on  the 
real  o-axis  and  represents  a  free  wave  Whose  phase  propagates  downstream  from 
the  excitation,  but  Whose  group  velocity  is  negative.  The  loci  4b  and  4c 
therefore  both  represent  waves  Whose  energy  propagates  towards  the  excitation 
point.  The  locus  labelled  Sc  terminates  on  the  real  o-axis  and  represents  a 
free  wave  which  propagates  downstream  from  the  excitation. 

The  difference  in  the  Character  of  the  response  at  frequencies  below 
0.128b  and  that  at  frequencies  above  0.168b,  together  with  the  form  of  the 
response  at  0.168b,  suggests  that  at  soam  frequency  between  these  values  there 
nay  be  a  branch  point  in  the  complex  w-plane,  as  a  result  of  poles  from 
opposite  half-spaces  of  the  complex  a-plane  coalescing.  The  form  of  the  root 
loci  shown  in  Figure  7,  show  that  such  a  pair  of  coalescing  poles  must  exist, 
and  hence  that  the  pipe  has  an  absolute  instability  which  exists  for  both  the 
impulse  excitation  and  the  time-harmonic  switch-on  excitation  whatever  the 
value  of  uq,  the  excitation  frequency.  Further  numerical  work,  not  presented 
here,  reveals  that  the  coalescing  poles  are  located  at 

Op  -  0.18-0.051  (4.1) 

<dp/2ir  -  0.1440.061  (4.2) 

Which  any  also  bo  obtainod  directly  from  equations  (3.4)  and  (3.5).  Tha 
component  of  the  pipe  response  due  to  the  integral  around  tha  branch  cut  in 
the  oonplmx  w-plane  grows  exponentially  with  tlxm,  and  at  larga  times  the 
response  is  approximately  at  the  frequency  wp  and  wavenumber  Op,  given  by 


equations  (4.1)  and  (4.2).  tes  frequency  ait  which  the  absolute  instability 
occurs,  aquation  (4.2),  doss  not  coincide  with  the  frequency,  obtained  from 
Figure  e,  at  teich  the  group  velocity  appears  to  be  Infinite. 

Piguree  8a  and  8b  ahov  the  real  and  Imaginary  parts,  respectively,  of 
the  wsvenueber- frequency  plot  for  the  undaspad  rubber  pipe,  obtained  free 
the  dispersion  relation  for  a  beam,  equation  (2.3).  Iha  branches  label  led 
1  and  2  are  purely  real,  and  branches  3  and  4  are  coupler  conjugate  roots, 
at  frequencies  between  fL  ( approximately  1.2BS)  and  fg  (approximately 
i.3Hs),  the  nueber  of  real  positive  roots  of  the  dispersion  relation  changes 
from  one  to  three.  At  frequencies  below  fj,  a  pair  of  coupler  conjugate 
roots  exist,  which  ooalaaoa  at  f-fL.  at  frequanclea  between  fj,  and  % 
all  four  roots  are  real.  At  £g  the  two  largest  real  roots  ooalaaoa,  and 
at  higher  frequanclea  this  pair  of  roots  are  ooaplex  conjugates. 

Figure  9  shows  a  root  locus  stability  plot  for  the  four  roots  ( labelled 
a,  b,  c  and  d)  of  the  und— part  rubber  pipe.  All  the  loci  label  led  1 
correspond  to  the  frequency  u^ZiHJ.TSBi.  Iha  loci  labelled  la  and  lb 
terminate  on  the  real  a-axis  and  represent  free  waves  dees  phase  propagates 
upstrexs  from  the  excitation  and  downstream  tram  infinity,  respectively,  the 
energy  of  both  of  these  waves  propagates  upstream  from  tea  excitation.  Hie 
locus  labelled  lc  crosses  the  real  e-axls  and  terminates  in  the  lower  half- 
plane.  Whilst  tee  locus  labelled  id  terminates  in  tee  upper  half-plane,  these 
represent  a  convective ly  unstable  wave  and  an  evanescent  wave,  reeectively, 
bote  propagating  downstream  from  tee  excitation,  the  loci  labelled  2 
correspond  to  tee  fkequency  «z/2ip-1.08b,  end  are  similar  to  those  at  0.79B8. 

the  loci  labelled  3  in  Figure  9  correspond  to  tee  frequency  wx/2*»l.2Mi 
the  loci  labelled  3a  and  3b  terminate  on  the  real  e-axls  and  represent  f ree¬ 
ve  ves  propagating  upstream  from  the  excitation,  tho  loci  labelled  3c  and  3d 
also  both  teradJiata  on  the  reel  a-axis  and  represent  free  waves  whose  phase 
propagates  downstream  from  the  excitation,  but  whose  energy  propagates 
upstream  from  Infinity  and  downstream  from  the  excitation,  respectively. 

The  loci  labelled  4  correspond  to  the  frequency  mz/2*»1.98b.  the  locus 
labelled  4a  terminates  on  tee  real  oMuds  and  therefore  represents  a  free- 
wave  propagating  upstream  from  the  excitation,  tea  locus  labelled  4b 
terminates  in  tea  lower  half-plane,  and  represents  a  wave  whose  amplitude 
decays  exponentially  with  distance  from  the  excitation  point,  tee  locus 
labelled  4c  terminates  on  the  real  o-exla  and  represents  a  free-wave 
propagating  downstream  from  the  excitation,  whilst  tho  locus  labelled  4d 
terminates  in  the  upper  half-plans  and  rsprssents  a  wave  whose  amplitude 
decays  exponentially  with  distance  from  the  excitation  point,  tee  waves 
represented  by  the  loci  4b  and  4d  have  positive  phase  speeds,  but  saro  group 
velocities,  tee  loci  labelled  5  correspond  to  the  frequency  ws/2n*1.79Bs, 
and  are  similar  to  those  at  l.SBs. 

The  only  coalescing  poles  tftiich  originate  In  opposite  half-planes 
ooalaaoa  at  a  real  frequency  between  1.298s  and  l.9Bs.  Zt  la  shown  In  ths 
Appendix  that  ooalesoenca  at  a  real  frequancy  corresponds  to  a  transient 
response,  and  the  rubber  pipe  therefore  has  no  absolute  Instability,  tee  pipe 
response  la  predominantly  at  the  excitation  frequancy,  and  tea  nature  of  tee 
response  depends  on  the  value  of  the  excitation  frequency,  if  the  excitation 
frequency  is  lsss  than  fj,  ths  pipe  Is  oonvsetlvsly  unstable,  if  the 
excitation  frequency  is  bstwssn  fj,  and  fg  ths  rsponss  consists  only  of 
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fg—  www,  and  if  the  —citation  frequency  ia  gr— tag  titan  f0  tha  raapon— 
la  again  atabla  and  oonsiata  only  of  fr—  —a,  an  eve— aoant  wave 
propagating  downatraaai  fr—  tha  excitation  and  a  —  who—  pha—  propagat— 
daw— tra—  fr—  infinity,  hut  who—  — gUlturta  dsoays  exponentially  with 
diatanoa  fr—  tha  axel tatien.  lha  valu—  of  ^  and  £p  any  ho  obtained 

aithar  toy  further  nuxerloal  work,  or  diroetly  fr—  aquatlona  (3.16)  and 
( 3.10).  For  tha  rubber  pipe  they  are  found  to  ha 

fx,  -  1.24Bb  (4.3) 

fg  »  1.2M1  (4.4) 

lha  m— rlral  roaulta  pro— ntad  in  Plgur—  7  and  9  ware  repeated  with 
pipe  dapping  included  in  tha  cm  leu  let  tone,  lha  nuna  rical  value  of  n»  the 
hyat erotic  lo—  factor,  w—  ehoaan  to  ha  0.2*  and  tha  raaulta  are  ahown  in 
Flgux—  10  and  U.  the  p— itiona  of  the  root  loci  for  the  at— 1  pipe,  shown 
in  figure  10,  are  only  a  lightly  altered  by  tha  addition  of  daaping.  In 
particular,  the  loci  oorraaponding  to  fr—  —a  propagating  upstreax  with 
nagativa  pha—  apaad  in  tha  undaapad  ca— ,  am  downwards,  i.a.  away  froa  the 
axis,  and  therefore  represent  evanescent  wav—  propagating  upetre— .  lha  loci 
corresponding  to  fr—  waves  propagating  upetrean  tram  the  excitation  with 
positive  pha—  epeed  in  tha  undaapad  ca— ,  nova  upwards,  i.a.  they  cro—  the 
real  at-axls,  and  therefore  rapes— nt  unstable  wav—  —  — plitudea  inarea— 
with  diatanoa  fr—  tha  excitation.  However,  it  is  clear  from  Figure  10  that  a 
locus  fr—  the  upper  half-pla—  nuat  c— lea—  with  a  locus  fr—  tha  lo— r 
half-plans  at  a  ooaplax  frequency  who—  real  part  is  bate— n  0.12—  and  0.1BHB 
and  that  tha  pips  is  therefore  absolutely  unstable.  Whilst  the  addition  of 
draping  to  the  pipe  — y  stablll—  or  deetabill—  individual  fr—  waves  title* 
propagate  at  the  excitation  frequency,  it  do—  not  raaovs  the  absolute 
instability.  Which  daainat—  the  solution  at  large  tin—. 

the  positions  of  the  root  loci  for  the  rubber  pipe.  Shown  In  Figure  11, 
are  again  only  slightly  altered  by  the  addition  of  daaping.  the  individual 
fr—  waves  propagating  at  the  excitation  frequency  — y  be  stabilised  or 
destabilised  in  the  sa—  way  —  for  the  at— 1  pipe.  However,  at  a  frequency 
of  approxlnately  1.2SHB  a  root  locus  fr—  the  upper  half-pla—  —alee— a  with 
a—  fr—  the  lower  half-pla—  at  a  coaplex  value  of  Wp,  indicating  an 
absolute  instability.  Whilst  the  addition  of  deeping  to  the  pipe  — y 

stablll—  or  deetabill—  the  individual  fr— . wav—  at  the  excitation  ftaquancy 

its  net  effect  is  to  d— tahili—  the  pipe  by  Introducing  an  absolute 
Instability  which  dosuLnat—  the  solution  at  large  tlx—  Irrespective  of  the 
excitation  frequency.  Thus  the  convective  instability  at  low  excitation 
frequenci—  is  converted  to  the  —re  serious  absolute  instability  which  occurs 
not  only  for  all  excitation  frequenci—  hut  also  for  iepul—  excitations . 

The  stability  characteristics  of  an  undaapad  pipe,  obtained  fr—  the  be— 
theory  nodal,  are  an— arieed  in  Figure  12.  If  8e/N>l  the  pipe  reap—  i« 
absolutely  unstable  for  both  the  lapul—  excitation  and  the  ti—  harxonic 
awltch-on.  if  aa/M«l  the  respon—  to  the  iapul—  excitation  is  stable,  being 
hounded  in  both  spa—  and  ti— .  The  xeepon—  to  the  ti— -harxonic  switch -on 
— y  he  oonvecti— ly  unstable  or  stable,  consisting  only  of  free  wav—  or  of  a 
aoeblna tlon  of  fr— -wav—  and  evanescent  waves,  depending  on  the  para— tare  of 
the  pipe,  tits  fluid  and  the  excitation  frequency. 


5.  C0MCL0DXM6  REMARKS 


•me  stability  characteristics  of  a  pips  con  vying  invlacld  fluid  hav 
boon  examirasd  using  a  approach  which  is  abla  to  distinguish  between 

an  absolute  instability,  in  which  the  amplitude  of  the  response  grows  with 
tine  everywhere,  and  a  convective  instabilty,  in  which  the  amplitude  of  the 
response  grows  as  the  disturbance  propagates  away  fron  the  excitation.  Two 
types  of  initial  disturbance  were  considered)  type  A,  whitfi  is  an  inpulse 
applied  at  t-O,  and  type  B,  which  is  a  tins  harnonlc  disturbance,  sin(UQt ), 
which  is  switched  on  at  t-O.  if  an  absolute  instability  is  present  the 
oscillatory  part  of  the  response  is  dominated  by  the  'pinch'  frequency.  Up, 
whilst  if  a  convective  instability  only  is  present  the  oscillatory  part  of 
the  response  is  predominantly  at  the  excitation  frequency,  uq. 

The  sain  findings  of  the  mathematical  and  numerical  «a>idc  on  the  beam 
theory  model  were: 

(i)  If  8m/M>l,  the  undamped  beam  has  an  absolute  instabilty  for  both 
type  A  and  type  B  disturbances.  One  mode  of  the  response  grows 
exponentially  with  time  everywhere,  and  eventually  dominates  the 
variety  of  other  wave  motions  which  may  also  be  present . 

(ii)  If  8m/IKl,  tbs  undamped  beam  response  to  type  A  disturbances  is 
transient  everywhere.  For  type  B  disturbances  the  form  of  the 
response  depends  on  the  excitation  frequency  as  follows: 

(a)  At  excitation  frequencies,  mq,  lew  than  a  small  valms,  *»{,, 
the  beam  is  oonvctivly  unstable.  Its  response  consists  of 
two  free  wave,  the  direction  of  whose  group  velocities  are 
upstream  from  the  excitation,  together  with  an  evanescent 
wav  and  a  oonvctivly  unstable  wav  whose  group  velocities 
are  downstream  from  the  excitation. 

(b)  At  excitation  frequencies,  uq,  in  a  range  to  wq,  the 

bean  la  stable.  Its  response  consists  of  four  free-mava,  two 
of  whose  group  vlocltlas  are  upstream  from  the  excitation  and 
the  other  two  of  which  are  downstream  from  the  excitation. 

(c)  At  excitation  frequencies ,  uq,  greater  than  wq,  the  beam  is 
also  stabla,  its  response  consisting  of  an  avanesoant  wav  and 
a  free  wav  whose  group  vloclty  is  upetream  from  the 
excitation,  and  an  evanescent  wav  and  a  free-wmv  whose 
group  vloclty  la  downstream  from  the  excitation. 

(ill)  The  addition  of  damping  cannot  rsemv  an  absolute  instability,  and 
it  introduces  one  fbr  both  type  A  and  type  B  excitations  if  one  is 
not  si  rawly  present,  the  e^onential  growth  factor  in  this  case 
being  proportional  to  the  loss -factor,  r\.  Because  the  losa-factor 
is  usually  email,  it  is  eiqpected  that  the  medlimi  term  behaviour  of 
the  beam  will  be  aa  described  in  ( i )  or  (ii)  abov,  but  the  long 
term  behaviour  must  be  absolutely  unstable. 
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anc*  will  b*  dsterained  by  non— linear  effects.  Thirdly,  the  results  presented 
here  are  for  infinite  pipes,  and  the  finite  length  of  practical  piping  irrstriM 
•ay  not  adait  the  long  wavelength  instabilities  predicted.  Finally,  the  de¬ 
stabilising  effect  of  daaping  aay  render  meaningless  the  analysis  of  the  un- 
daaped  case  because  some  daaping  is  always  present  in  practice.  The  best  way 
to  test  the  range  of  applicability  of  the  theory  is  to  conduct  experimental 
tests. 


E.  h.  Skelton  (HSO) 

Manuscript  completed  January  1984. 
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D  X  X 


A  causal  Approach  to  the  stability  Analysis  of  Mnsar  fyet  — 


In  tha  ataady  atata  approach  to  atabllity  it  i s  assumed  that  aa 
excitation  of  the  fora  Be[Fo«xp( -iugt )]  exists  for  all  time  t.  Zn 
practice,  hcwater,  tha  exaltation  will  hase  bean  switched  on  at  same  ti me, 
t-O,  say.  throe  possibilities  for  the  aapUtude  of  tha  system  response  than 
arise  i  (1)  it  nay  settle  down  everywhere  to  finite  values  which  are  bounded 
in  both  space  and  tiro  and  which  nay  or  nay  not  be  equal  to  the  steady  state 
values i  (li)  it  nay  settle  down  mrywhere  to  steady  values  Which  increase 
exponentially  with  distance  from  the  point  of  excitation;  or  (111)  It  nay 
grow  exponentially  with  tlxn  at  every  point  of  the  systsn.  Ms  letter  [11]  and 
others  have  tewed  ease  (11)  a  convective  instability,  and  case  (ill)  an 
absolute  instability.  The  stability  analysis  proosade  as  follows. 

It  is  convenient  to  represent  the  system  response,  x^a.t ),  and  the 
excitation  r,(s,t)  as  the  real  parts  of  coupler  functions  W(s,t )  and 
P(x,t),  respectively,  via.. 


These  fUnctlone 


wr(x,t)-»ats(a,t)] 

Pr(x,t)-iteCP(>.t)] 


(hi) 


than  represented  as  the  Taplars  rowrter  transform  pair 


ob+Io  m 

*s.t) 

2  f  f 

i(o,e) 

-  (1/de  ; 

P(a.t) 

J  J 

?(a,e) 

-eWfl  -e 

exp(  laa-lwt  )dadM 


ii(e,e) 

P(«,e) 


» 


I 


»(s,t) 


exp(  -last tut  )dadt 


|F(s,t) 


(A2) 


(A3) 


whsro  a,  a  positive  constant,  is  chosen  to  be  such  that  the  line  in  the 
complex  e-plane  rm[e]-c  lies  above  any  singularities  of 


S(s,e)-(l/2*) 


|  S(a,w)exp(  lam)da 


(AA) 


This  ensuroe  that  causality  la  satisfied,  l.e.  the  response  of  the  systaai  is 
aero  for  t<0. 


d 


/  v  m 


(AS) 


The  Laplace-Fourier  trtMfoni  of  the  linear  differential  equation 
satisfied  by  the  response  W(z,t)  My  bo  wrttton  as 

D(a,w)W(a,<*)-F(a,w) 

mwzt  D(a,u)  is  tbs  dispersion  eolation  of  tho  systea.  substituting 
oquatlon  (AS)  into  aquation  (A2)  givoo  tho  syotoa  rooponoo  as 

•4-iS  CB 

0(s,t)-(V4ir2)  j  |  exp(las-ist)d«dw  (AS) 

-OBfio  -ao 

Tho  first  stop  in  tho  evaluation  of  the  response  W(  z,t ),  is  tho 
evaluation  of  tho  Fourier  (w.r.t.  a)  integral  of  aquation  (AS ),  which 
requires  intonation  about  the  behaviour  of  both  D(a,u)  and  F(a,w).  The 
dispersion  relation  D(a,u)  is,  in  general,  a  sulti-valued  function  >*iiCh 
ny  be  sate  single-valued  on  the  real  a-axis  by  imposing  any  radiation 
condition.  In  any  probleM  involving  fluid-structure  interactions,  for 
exaspla,  the  dispersion  relation  involves  the  radical  vC(e-o0)2/c*-a2] 
which  My  be  node  single  valued  by  iapoaing  the  radiation  condition  that 
Ia(y((u-oO)2/c2-a2]]  Should  be  non-negative  on  the  real  a-axis.  The 
(single-valued)  dispersion  relation  on  the  real  a-axis  My  then  be 
analytically  continued  into  the  complex  a-plane,  where,  in  order  to  ensure 
that  the  continuation  is  single-valued  everywhere,  branch  cuts  oust  be 
introduced.  For  sufficiently  large  a  the  branch  points  will  be  off  the 
real  a-axis.  The  aseoclated  branch  cuts  aust  not  cross  the  real  oi-axie 
and  aust  therefore  go  to  infinity  in  their  respective  half-planes.  With  the 
exception  of  this  restriction,  the  Choice  of  branch  cuts  is  arbitrary.  The 
behaviour  of  F(a,w)  depends  only  on  the  type  of  forcing  excitation  applied 
to  the  systea.  Two  ooaann  types  of  excitation  axe  an  iapulslve  excitation 

F(z,t)-F08(x)6(t)  (A7) 

whose  tranefOxx  is  given  by  equation  (A3)  as 

F(a,«)-F0  (AS) 

and  a  tiae-hamonic  excitation  switched  on  at  tr-O,  vis.. 


F(z,t)- 


o 


t<0 


F0  8(z)exp(-lw0t )  t»0 


(A®) 


whoee  transform  ie 


F(«,»)-iF0/(<i»-«0) 


(A10 


The  Fourier- integral  of  equation  (AS)  aay  be  evaluated  by  suitably 
closing  the  contour  and  applying  Cauchy's  theorem.  Figure  2  illustrates 
the  case  tdien  z>0,  requiring  that  the  contour  is  closed  in  the  upper  half¬ 
plane.  for  s<0  the  contour  aust  be  closed  in  the  lower  half-plane.  This 
convention  ensures  that,  as  Za(a]  tends  to  plus  or  minus  infinity. 
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respectively,  the  term  exp(iaz)  in  the  integrand  remains  finite. 

in  Pigure  2,  let  ( j-l,Na)  represent  all  the  poles  of  the  integrand 

in  the  upper  half-plane,  and  aj  ( j«l,Na)  represent  all  the  poles  in  the 
lower  half-plane.  If  a  is  sufficiently  large  then  all  these  poles  are  well 
separated  from  the  real  a- axis .  The  contours  and  C2  are  circular 

area  of  large  radius  and ,  provided  that  all  these  singularities  of  the 
integrand  in  the  upper  half -plane  lie  inside  the  contours,  the  contributions 
to  the  integral  from  integrating  along  and  C2  tend  to  zero  as  the 

radius  tends  to  infinity.  The  integrals  along  B1,  b2,  and  b3  are  branch 
cut  integrals  and  will  be  discussed  later.  Similar  results  are  true  for  the 
contour  for  z<0  in  the  lower  half-plane.  The  integral  may  be  generally 
expressed,  as  the  sum  of  residues  and  branch  cut  integrals,  in  the  form 


W(z,w)  -  *i ^  P( a* , «i) )exp(  ia*z )/X)  ( a* , w )  +  B  (z,w)  (All) 

J  J  J 

Where  the  +  sign  corresponds  to  positive  values  of  z  and  the  -  sign  to 
negative  values  of  z.  Ba(z,u)  is  the  contribution  from  the  integral  around 
the  branch  cuts  in  the  complex  a-plane,  and 

i 

D  (a,u)-dD(a,u)/da 

The  system  response  W(z,t)  is  then  found  by  evaluating  the  Laplace- 
integral  of  equation  (A6),  using  the  result  of  the  Pourier- integration  given 
by  equation  (All).  Por  t<0,  the  contour  must  be  closed  as  shown  in  Pigure  3a 
in  order  to  satisfy  the  causality  requirement.  Por  t>0,  the  contour  may  be 
closed  by  using  the  standard  lap lace  contour  also  shown  in  Pigure  3a.  Any 
numerical  scheme  for  evaluating  this  Laplace  integral  would  become  numerically 
unstable  for  large  positive  values  of  t  because  of  the  term,  exp(ot),  in  the 
integrand  Which  tends  to  infinity  as  t  tends  to  infinity  for  positive  values 
of  a.  it  is  therefore  necessary  to  deform  the  contour  down  to  the  real 
u-axis  by  letting  cr  tend  to  zero,  not  only  to  allow  a  numerical  evaluation 
of  the  integral,  but  also  in  order  to  identify  the  contributions  to  the 
integral  from  branch  cuts  and  poles  in  the  complex  u-plane.  The  resulting 
contour,  shown  in  Pigure  3b,  lies  just  below  the  real  w-axis,  but  is  indented 
around  any  poles  and  branch  points  in  the  upper  half-plane  or  on  the  real 
u-axis.  If  un  (n**o,Nu)  represent  all  the  poles  of  the  integrand,  then  the 
system  response,  for  t>0,  may  be  expressed  using  Cauchy's  theorem,  as  the  sum 
of  residues  and  branch  cut  integrals,  in  the  form 

N  N 

u)  a 

W(z,t)  -  £{-iBa(z,wn)±£  R<a*,«n)exp(ia*z))exp(-iwnt)  +  Bw(z,t)  (A12) 
n-O  5*1 

where  R(a*,u  )  and  B  (z,u  )  are  the  residue  contributions,  defined  as 
1  n  an 

R<a*,w  )  -  lim{(  w-w  )P( q* ,u)/d'  ( a* , w ) }  (A13) 

*  n  u-u  n  1  J 
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B  ( *,<*> )-li*{(<*-ti  )B  ( *,«)} 
a  n  n  a 


(A14) 


and  Bu( x , t )  is  the  contribution  from  the  integral  around  ths  branch  cuts 
in  ths  complex  o-plane. 

In  order  to  deform  ths  contour  in  ths  u-plane  as  described  above,  ths 
continuation  of  z,u)  into  ths  region  In[u]<a  oust  bs  analytic,  as 
ths  Laplace  contour  is  Moved  down  to  ths  real  u-axi a,  ths  poles  enclosed  by 
ths  Fourier  contours  also  Move  in  ths  o-plans  and  it  nay  be  necessary  to 
defora  the  Fourier  contour*  so  that  they  always  include  the  sane  poles  of 
the  response  transfora,  w(a,w). 

in  order  to  ascertain  whether  or  not  any  physical  significance  can  be 
attached  to  each  tern  of  equation  (A12)  for  the  response  of  the  systea,  it 
is  necessary  to  use  the  root  locus  technique  as  described  by  Mb  letter  [11] 
and  Atkins  [14].  Firstly,  consider  the  contribution  fzoa  a  staple  pole  on 
the  real  u-axis  at  wwq,  which  arises  froa  a  tine— haraonic  excitation  of 
the  fora  F0exp(  -i^t ) .  In  the  root  locus  technique  the  loci  of  roots  a| 
(which  are  functions  of  w)  are  plotted  for  values  of  u-Wy+io  for  particular 
values  of  Up  as  a  varies  froa  4m  to  o.  Provided  that  these  loci  do  not 
coalesce  (i.e.  intersect  at  the  ease  value  of  Wf  and  a)  for  any  real 
then  there  are  twelve  basic  possibilities  for  each  locus  which  are  illustrated 
in  Figure  4a.  Ihe  shaded  circles  in  this  Figure  denote  points  where  the  loci 
terminate,  corresponding  to  o*0. 

the  fora  of  the  systea  response  aay  be  detezained  froa  the  loci  for  which 
er*Re[un].  In  cases  (1)  to  (6),  because  each  of  the  loci  coaes  froa  above  the 
real  o-axls  the  Fourier  contour  must  be  indented  where  necessary  to  include 
the  pole  in  the  upper  contour.  The  response  due  to  eadi  of  theee  poles  is 
non— rero  only  for  s>0,  in  which  case  it  is  given  by 

M(s,t)-F0erp{(ias-a1)s}eicp(-iunt)  (A15) 

where  a +ia.  is  the  location  of  the  pole  when  o*la[u  ]-o. 


In  cases  ( 1 )-( 3 ) ,  where  as>0,  the  phase  disturbance  propagates  in  the 
positive  s-dlrectlon.  If  a^>0,  as  in  case  (1),  the  wave  decays  exponentially 
with  distance  and  it  is  called  an  evanescent  wave,  if  o^-O,  as  in  case  (2), 
the  wave  propagates  without  decay.  In  case  (3),  where  a^<0,  the  wave 
increases  exponentially  with  a  and  is  tamed  ' convectively  unstable'.  In 
cases  ( 4 )-< 6 ) , where  «s<0,  the  phase  disturbance  propagates  in  the  negative 
s-directlon  (i.e.  towards  the  excitation  ).  In  case  (4),  where  a^>0,  as  the 
phase  propagates  the  wave  amplitude  increases  exponentially  with  distance. 

In  case  (5)  o^-o  and  the  wave  propagates  without  decay.  If  a^<0,  as  in 
case  (6),  the  wave  decays  exponentially  with  distance. 

cases  ( 7 )— ( 12 )  are  similar  to  cases  ( 1 )-< 6 ),  except  that  the  loci  coaw 
from  below  the  real  a-axis  and  the  Fourier  contours  must  be  indented  where 
necessary  to  ensure  that,  in  each  case,  these  poles  are  included  in  the  lower 
contour.  She  contribution  to  the  response  from  each  of  these  poles  is 
therefore  non-sero  only  for  s<0. 


The  contributions  to  ths  displacement  fro®  simple  poles  in  the  complex 
u-plane  represent  propagating,  spatially  amplifying  or  spatially  decaying 
waves  at  the  excitation  frequency,  This  description  of  the  propagation 

directions  gives  only  the  phase  speeds  of  the  individual  waves.  A  more 
realistic  quantity  physically  is  the  group  velocity,  du/da,  the  velocity  at 
which  energy  propagates,  this  may  be  obtained  from  the  direction  of  the 
loci  at  <T-0,  viz., 

dw/acp^icV^a^+ic^J-Cfdfl^/dsHifda  r/9a)}(d<r/da)  (A16) 

where  s  is  the  arc  length  along  the  locus,  measured  from  o» 0. 

In  particular,  the  group  velocity  is  positive  if  the  initial  direction  of 
the  locus  in  the  complex  ot-plane  is  upward,  and  the  group  velocity  is 
negative  if  the  initial  direction  of  the  locus  is  downward. 

In  order  to  consider  the  contribution,  Bu(z,t),  to  the  response  from 
the  integrals  around  the  branch  cuts  in  the  complex  u-plane  possible 
singularities  in  the  complex  u-plane  must  be  identified.  Equation  (All) 
indicates  that  singularities  in  the  u-plane  may  be  expected  if,  for  some  u 


D  («j(u),u)-0 


(A17) 


where 


«2<«) 


is  a  root  of  the  dispersion  relation 


D(a^(u),u)-0 


(A18) 


If  a  solution  to  the  simultaneous  equations  (A17)  and  (A18)  exists,  with 
Im[u]>o,  it  indicates  that  equation  (A18)  has  a  repeated  root  aj  at  that 
value  of  u,  which  in  the  context  of  the  root-locus  method  means  that  two 
root  loci  coalesce  in  the  complex  a-plane.  There  are  two  distinct 
possibilities  here  as  illustrated  in  Figure  4b.  If  both  loci  have  come  from 
either  above  or  below  the  real  a-axis,  as  in  case  (l),  then  the  Fourier 
contour  in  the  complex  a-plane  may  always  be  deformed  to  include  both  poles, 
and  at  the  coalescence  frequency  the  residue  may  be  calculated  as  for  a  second 
order  pole  at  ot-a j  ( u ) :  this  constitutes  a  removable  singularity  in  the 
complex  u-plane.  Thus,  as  pointed  out  by  Helcher  [11],  the  coalescence  of 
poles  illustrated  in  case  (1)  does  not  indicate  a  true  singularity  in  the 
complex  u-plane.  In  case  (2),  however,  a  locus  originating  above  the  real 
a-axis  coalesces  with  one  originating  from  below.  The  Fourier  contours  in  the 
complex  a-plane  can  no  longer  be  distorted  to  include  only  the  poles 
originating  in  the  upper  (or  lower)  half-plane,  and  they  become  'pinched'  at 
the  (complex)  frequency.  Up,  at  which  coalescence  occurs.  This  confirms  the 
presence  of  a  singularity  in  the  complex  u-plane  at  u-%>p,  which  may  be 
identified  as  a  branch  point  by  considering  the  Taylor  series  expansion  of  the 
dispersion  relation  about  u-up,  a-Op,  viz.. 


D(a,w)“(w-w  )dD/aw+( l/2)(a-a  )282D/8a2+(a-a  )(u-u  )d2D/dadw 
P  P  P  P 


+(  1/2  )(  u-u  )282D/8u2+ .  (  A19 ) 

P 

where  all  the  derivatives  are  to  be  evaluated  at  a-ap,  u-wp.  By  considering 
the  order  of  each  of  the  terms  in  equation  ( A19 )  as  u-up,  Helcher  [11]  has 
shown  that  in  a  neighbourhood  of  (Op,up) 
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(MO) 


( ct-a  >-±[2(  a>-a>  )(  «D/9u)/(  *2D/«o2  )  ] 1/2 

P  P 

The  residue  contribution  to  W(  z,a>)  in  the  neighbourhood  of  (ap,Up),  la  giv«m 
by  aquation  (All)  aa 


*iF(a,u)erp(ias)/D>(a,<ri)  (A21) 

Diffarantiation  of  tha  Taylor  expansion,  aquation  (A19),  with  raapact  to  a 
allowa  tha  raaidua  contribution,  aquation  (A21),  to  ba  written  to  firat 
order  in  (c*-ap)  aa 


*iF(a,u)exp(las)/((«-a  )92D/«a2)  (A22) 

P 

which  may  be  expressed  in  tens  of  m,  uaing  equation  (A20),  a a 

*if(a,M)eaqKia*)/t2(«-“p>(a2D/a«2)(  9D/9m)]l/2  (A23) 

dmeonstrating  that  at  is,  in  fact,  a  branch  point  of  W(  * ,  o> ) . 

P 

In  the  context  of  the  root-locus  method,  the  presence  of  the  branch 
cut  may  be  demonstrated  by  considering  tha  possible  behaviour  of  neighbouring 
root- loci  in  the  complex  <x-plane  for  ot_  and  «+,  where  ut_<Re[atp)  and 
«+>»•[>,,],  as  illustrated  in  Figure  5a  and  5b.  In  both  cases,  jit  may  be 
Shown  mat  the  residue  contributions  to  the  Fourier  integral  W(s,u)  are 
discontinuous  at  Re(u]-Ra[Op],  indicating  that  the  branch  cut  in  the  complex 
M-plane  runs  from  otp  to  Op-im.  this  Choice  of  branch  cut  is  imposed  by 
Idle  causality  requirement  which  determines  the  poles  to  be  included  in  each 
Fourier  contour. 

The  approximate  value  of  the  contribution  to  the  response  from  the 
integration  around  the  branch  cut  in  the  complex  u-plane  may  be  calculated 
by  miking  the  assumption  that  F(a,u)  may  be  approximated  by  F(ap,up). 

This  approximation  is  exact  for  an  impulsive  excitation,  and  the  branch  cut 
integral  then  becasms 


Bja.t)  - 


±iT(  «p,  )«*P(  l<*p* ) 

tr(2(  92D/9a2 )(  90/ )  ] 1/2 


which  can  be  evaluated  analytically  as 
B^(*.t)  - 


exp<  -io)t ) 


,1/2 


dot 


u-im 

P 


±U  1-i  )F(  ap,  up  )exp(  jgps-iwpt ) 
2[ir( »2D/«a2 X 90/ 9m ) ] 1/2 t 1/2 


(A24) 


(A25) 


Aquation  (A25)  Shows  that  if  the  imaginary  part  of  up  is  positive 


UNLIMITED 

the  contribution  to  the  response  from  the  branch  cut  integral  grows 
exponentially  with  tins  at  all  locations,  a.  this  is  termed  an  abaoluta 
instability,  and  occurs  whatsvar  tha  valua  o£  the  excitation  frequency, 
and  in  fact  it  also  occurs  for  an  impulsive  excitation.  When  a  aystaai  is 
absolutely  unstable,  tha  systaa  response  is  always  dominated  by  tha  ’pinch' 
frequency  and  wavenusber,  Up  and  Op  respectively,  after  a  sufficiently 
long  tine  interval.  If  Up  is  real  then  the  factor  t-1/2  in  equation 
( A25 )  ensures  that  the  response  is  transient. 

It  is  now  necessary  to  discuss  the  branch  cut  integrals  in  the  cosplex 
a— plane .  in  general,  provided  that  the  brand)  cuts  are  well  away  froai  any 
singularities  of  the  integrand,  the  integral  around  the  circular  contour  B 
is  zero,  and  the  integrals  along  the  contours  Bj  and  B3  contribute  only 
to  the  near  fields  because  the  integrands  are  bounded.  Crighton  [3]  has 
pointed  out  that  only  those  poles  which  exist  for  all  suitable  choices  of  the 
branch  cuts  nay  have  any  physical  significance  at  large  values  of  z.  The 
regaining  poles  contribute  only  to  the  near  field  and  their  contributions 
must  always  be  combined  with  those  frosi  the  branch  cut  integrals.  The  precise 
form  of  the  branch  cut  integrals  depends  on  the  choice  of  branch  cuts  and  on 
the  form  of  the  dispersion  relation. 

As  already  mentioned,  if  a  is  sufficiently  large  all  the  poles  of  the 
integrand  of  equation  (A6)  are  well  separated  from  the  real  o-axis.  A 
consequence  of  this  is  that  the  location  of  the  poles  depends  on  a  particular 
choice  of  branch  cuts,  and  that  no  pole  exists  for  all  possible  Choices  of 
branch  cuts.  It  also  means  that,  as  a  varies,  a  root  locus  corresponding 
to  ths  movement  of  a  particular  pole  may  cross  a  branch  cut  and  ’disappear' , 
or  vice-versa.  At  a-O,  the  poles  corresponding  to  evanescent  waves  also 
depend  on  the  Choice  of  branch  cuts.  Ths  contribution  to  the  response  from 
these  poles  can  therefore  have  no  physical  significance  in  isolation,  but 
must  be  coablned  with  the  branch  cut  integrals.  A  free-wave  pole  (i.e.  one 
Which  lies  on  the  real  a-axLs  for  <r-0 )  always  exists  for  all  possible 
choices  of  branch  cut  and  therefore  has  physical  significance.  A  pole 
corresponding  to  a  convective  instability  crosses  the  real  o-axis  at  some 
positive  value  of  a  and  therefore  always  exists  in  a  neighbourhood  of  the 
axis.  Furthermore,  according  to  Atkins  [14],  if  coalescing  poles  in  the 
a-plane  corresponding  to  a  singularity  in  the  upper-half  w-plane  exist,  they 
do  so  for  all  possible  Choices  of  branch  cut. 

In  the  absence  of  branch  cuts  in  the  complex  a-plane,  however,  the 
location  of  all  the  complex  poles  in  the  a-plane  may  be  uniquely  determined, 
and  their  contributions  to  the  response  may  be  identified  as  propagating, 
evanescent  or  oonvectively  amplifying  waves. 

Vo  si  swim  lie,  the  root  locus  technique  may  be  used  to  determine  ths 
stability  character 1st ics  of  a  linear  system.  Firstly,  it  is  necessary 
i  to  plot  the  root  loci  for  a  range  of  frequencies  in  order  to  determine 

J  whether  or  not  there  is  a  branch  point  in  the  complex  u-plane  as  a  result 

of  coalescing  poles  Whose  loci  in  the  a-plane  cosw  from  opposite  half-planes . 
If  such  a  branch  point  exists  the  system  Is  absolutely  unstable  and  the 
e  response,  which  is  predominantly  at  the  branch  point  frequency  Wp  grows 

I  exponentially  with  time  at  all  locations.  If  no  such  branch  point  exists, 

the  response  of  the  system  may  be  determined  by  the  root  loci  in  the  a-plane 

i 
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which  correspond  to  the  excitation  fnqumelM.  .  nodes  of  the  response  nay 
bs  identified  as  non— decaying  travelling  waves  or  owveetlvely  aaplifying 
waves,  and,  in  cases  whsre  there  are  no  branch  points  in  ths  a-plane, 
evanescent  waves  nay  also  be  Identified,  in  practice  Instabilities  are 
usually  ilmp ail  by  soaa  non-linear  mchanisn,  or  other  factors  not  taken  Into 
account  in  the  current  theory. 

A  necessary,  but  not  sufficient,  condition  for  Instability  (absolute 
or  convective)  is  that  a  locus  In  ths  a-plane  crosses  ths  real  at-axis. 

This  laplies  that  for  soaa  real  wavenuntoer  a,  ths  dispersion  relation  oust 
have  a  root  u  with  a  positive  imaginary  part,  nils  is  ths  criterion  which 
is  usually  used  to  determine  instability  using  the  steady  state  approach. 


